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Inlroduction 


The curve which passes through two specified points with specified orientation while minimizing 


S 


= /« 


ds, 


where k is the curvature and s the arc-distance, has a number of interesting applications. 

In a thin beam, curvature at a point is proportional to the bending moment [pg. 80,1], The total elastic energy 
stored in a thin beam is therefore proportional to the integral of the square of the curvature [pg. 163,1], The 
shape taken on by a thin beam is the one which minimizes its internal strain energy. This is why we call the 
curve sought here die minimum energy curve. A thin metal or wooden strip used by a drafts(wo)man to 
smoothly connect a number of points is called a spline [pg. 156, 2], Such splines are used in creating lofted 
surfaces from plane parallel cross-sections of ship hulls and aircraft fuselages [pg. 228, 2]. The shape of a 
spline constrained to pass through two specified points with specified orientation is what we are after here. 


In computer graphics and computer aided design there is a search for curves which are particularly "smooth" 
[pg. 156, 2; pg. 49, 3; pg. 119, 4; pg. 309, 5; pg. 43, 6]. One measure of smoothness is die inverse of the 
integral given above. Typically, cubic polynomial approximations are used instead of the optimal curve 
[pg. 162,2; pg. 66, 3; pg. 129, 4; pg. 315, 5]. (Unfortunately, these approximations arc called splines too). 

It has been suggested that the human visual system uses a curve of low energy when completing a contour. 
Ullman [pg. 1, 7] proposes that a subjective contour consists of two circular arcs tangent at dieir point of 
contact. Out of the one-parameter family of solutions of dais form he picks the one which minimizes the 
integral of the square of curvature [pg. 2, 7], He notes diat die curve so constructed may have near-minimal 
energy. (This does not mean that it necessarily lies close to die curve of minimum energy, as we shall see.) 
Brady et al used cubic polynomial approximations instead [8]. 


Preview 

We will first consider a special case. Here die curve must pass through the points (-1,0) and ( + 1,0) in the 
jcrplanc with vertical orientation at both points. We first determine the optimum curve as the limit of a series 
of approximations. This approach helps suggest algorithms for computing approximations to the ideal curve. 
Later we solve the variational problem directly. We dicn discuss how this curve can be translated, rotated, 
and scaled to produce a four parameter family of curves which contains the general solution. 


A Semicircle 

One curve which connects the two points (-1,0) and (+1,0), and has the desired orientation at these points, is 
a semicircle of radius one with the center at the origin (see Figure 1). The curvature equals one at all points 
and so the relevant integral for the section in the right hand quadrant becomes 



1 ds = m/2. 


The arc length happens to have the same value, 


5’- 



ds = m/2, 


while the maximum height of the curve above the jc-axis is 


% = 1. 

Can we do better, that is, find a curve with a smaller value of 8? 


Two Circular Arcs 

We try a combination of two circular arcs for the portion of the curve in the first quadrant, as shown in 
Figure 2. (The other half of the curve is obtained by reflection about the y-axis.) Let the first arc have radius 
R and angular extend 6 , while the remaining portion has radius r and angular extend (tt/2-0). Note that the 
parameters r, R, and 6 are not independent, since one can obtain from the diagram 

(R -r) cos 6 = (R-l). 

The arc length of the right-hand portion of the curve becomes 

= 6 R + (m/2 -0) r, 

while the energy is 8 = 4 + . 

R r 

Minimizing 8, subject to the given constraint, using the method of Lagrangian multipliers, leads to the set of 
equations 


1 

R 


0 


R 


2 + ^ : ^ (l-sec 0 ) = 0, 


l + fe/ L:<7)(* —r) tang = 0) 


(R-r) cosfl - (R-l) = 0. 


ITte second of these can be simplified into 


r - R (tt/2 -0) tan 6 , 

which, when applied to the first equation, yields 


Solving this numerically we obtain 
We can also show that 

so that 


0 (w/2 -0) (1 + sec 0) — 1. 

0 = 0.412868765... 

: _ 1 _ 

[(1 - cos 0) + (ir/2 -0) sin 0] ’ 

R = 1.8227161... 


and 


r = 0.92452847... 



Finally, since here 
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we see tliat 


and 


while 


= r + (R -- r) sin 6 


3G = 1.2849161... 


g = 1.4789649... ~ 0.94153834*(w/2), 
S= 1.8230795... ~ 1.1606084*(ir/2). 


'1’he energy in this curve is only 94.15% of that in the semicircle, so we can do better. 


The Best Ellipse 

The two-arc solution suggests that the optimum curve is elongated and has radius of curvature smaller than 
one at its peak, and larger than one near the .v-axis. Is it an ellipse? 1 he equation of the ellipse [pg. 411,13, 

pg. 72.14] shown in Figure 3 is 


or in parametric form, 


(. x/a) z + {} ! /b) z - 1, 
x - a cos /, y - b sin /. 


The eccentricity e is defined by the equation 


> 2 = 1 — ( a/b) z . 


The arc length can be found as follows [pg. 26,9], 


if = J (is = f* /Z -y/ x 2 + / dt 


where 


Now, 


or. 


dx . • dy 


dt ‘ 


x 2 + j 2 = < 2 2 sin 2 f + b z cos Z U 


£ + y 2 = b z [1 - e z sin 2 /]. 


So 


J = b r n y[T-e z sin 2 / dt = b E(e), 


where E{e) is die complete elliptic integral of the second kind [pg. 16,9; pg. 833,10; pg. 904,11, pg. 589,12]. 
(The elliptic integrals got their name from the fact that they first appeared in the mensuration of the ellipse.) 

Finally, for a= 1, 


= m l ■/ 


1-e 2 . 


1. Note that E denotes two things in this paper. When it has one argument,as here, it signifies the complete 
elliptic integral of the second kind. 



The curvature can be found as follows [pg. 553,13; pg. 22,14]: 


K 


• « • • • • 

y * ~ 


ab 


[i 2 + ; ,2 ] 3/2 [a 2 sin 2 / + b 2 cos 2 /] 3/2 ’ 


Now 

8 = Jk 2 ds = f k 2 [i 2 + y 2 ] i/ 2 dt, 

So 

- a 2 p /2 1 , 

b 3 Jo [1 - e 2 sin 2 /] 522 


r l’he definite integral can be shown [using pg, 165,11] to be equal to 

[2(e 2 + 1) F.(e) - e' 2 K(e)}/0 A 
where e 2 + e' 2 = 1, 

and K(e) is the complete elliptic integral of the first kind [pg. 16, 9; pg. 834, 10; pg. 904, 11; pg. 589, 12]. So 
finally when a= 1, 

8 = [2 (2 -e 2 ) E(e) - (1 -e 2 ) K(e)) / [3 1-e 2 ]. 

When we set e=0 we obtain & = tt/2, as we should, since the curve in this case is just a semicircle. 

In order to find the best possible ellipse we need to differentiate the expression above with respect to c. Here 
we need the following derivatives [pg. 21, 9; pg. 907,11]: 


dm - E(k) - K(k) 

dk ~ k 

dim E{k ) - k' 2 Ml 

dk k k' 2 


where k 2 + k' 2 - 1. 

The eccentricity of the optimum ellipse satisfies 

E(e) [4e A - 5e 2 + 3] = K(e) [2c 4 - 5e 2 + 3]. 
Solving this equation numerically leads to 

e = 0.6530018... 

% = b/a = 1 / 1-c 2 = 1.3203823... 

8 = 1.4674751... ^ 0.93422368 * (w/2). 


with 


■ 6 ■ 

y = 1.8311202... « 1.1657273 *(v/2). 

The maximum and minimum radii of curvature are 

r = b z /a = 1.7434096... and r . = a 2 /6 = 0.75735636... 

max min 

This curve has an energy which is only 93.42% of that of the semicircle. We have found a curve which has a 
smaller value of 8 than our two-arc solution. 

Can we do better still? 


Multi-arc Approximation 

Consider a smooth curve constructed out of n circular arcs (see Figure 4). Let the radius of curvature of the 
piece turning through the angle from a. to « i+1 be r.. We note that the total arc length, 3", and the integral 

of the square of the curvature, 8 ', are given by 


n-1 m a 

*'=2 «i> r i = Vl«n + 2 (Vl- r iK’ 

i *0 


n-1 

2 

i a l 


i «0 1 


K 

r n-l i = l Lr i 'i-l J 


where a Q = 0 and a„ = w/2. We also have to compute the width, V', and height, % ', 


n 


n-1 


n-1 


r = 

2 

r. (cos - cos a i+1 ) : 

= >0 cos « 0 - 2 (^i-i - 


i *0 

i = l 


n-1 


n-1 

%' = 

2 

r. (sin a i+1 - sin a .) = 

Vl sina n +2 


i*0 

i^l 


To solve our original problem we need to scale whatever curve we obtain so that its overall width equals 2 
instead of 2WThe scaled values are as follows: 

% = % vr, 

8 = g'lT, 

y = r/w. 

Note that tire integral of the square of the curvature is decreased when we make the curve larger without 
changing its shape. 


t 
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Optimum Multi-arc Approximation 

Our task now is clear: We have to minimize 6 ' t' by suitable choices of the parameters r. (for /= 0,1... n -1) 
and a i (for /= 1,2 ... n -1). Actually, we can pick one of the radii arbitrarily, r for example, since the whole 

curve will simply be scaled accordingly. I he minimization looks difficult at first when one considers the 

complexity of the product 8 ' If' and its derivatives. It appears necessary to resort to numerical techniques to 
solve for the 7(n -1) parameters. 


Fortunately this is not the case, for if 


^(8'r) = o, 

then, by the rule for the differentiation of a product. 


9Tf' .98' 
9x ' 9x 


If 78', 


for arbitrary x 0- e. r. and a.). Since the right-hand side is independent of x, it must equal a (positive) 
constant, c 2 say. Thus we find that 


91f' ,98' 2 


/ 


3a i ' 3a. 


c for / = 1,2 ... n -1, 



air. 88' 2 p , 

877/ jr =c for, = 0.1... «-l. 


We need the following derivatives now, 


98' 

9a i 

98' 
9 r. 


91f' 


ri_ 

. i ] 


r i-i J 


(a. +1 - a.) 

- ~2 -for i = 0,1... n -1, 


foT = O',-! “ r i )sin a i for i = 1,2 ...«-1, 


91f' 


= (cos a i - cos a i+1 ) for / = 0,1... n -1. 


Using these derivatives in the equations above we obtain 

r. sin a. = c 2 for / = 1,2... n -1, 


cos a. . - cos a, , 

_ jZ -i±l-1 = c 2 f or i - qj ... „ -i 

1 «i + i ~ «i 
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It is easy to verify that in the ease that n- 2, we obtain the same equations as before, provided we introduce 
the additional constraint W'=l or 

r o “ < r o-- V cos af i 

Note that, if c is known, a simple procedure will give us all of the parameters. Let r Q = 1, say, then the second 
equation can be used to find a.. (This non-linear equation has to be solved numerically.) The first equation 
then allows one to solve for r Knowing r r the second equation allows one to find a v and so on. If the 
value of c is correct, the process will terminate with a n = ir/2. The correct solution can be found by 
searching for the appropriate value of c. This is very much simpler than a direct search on the 2 (n -1) 
parameters. 


Sonic Helpful Relationships 

A number of interesting observations can be made now about the multi-arc solution. First of all, the "energy" 
68 ' in an individual arc is directly proportional to the projection 5W' of this arc on the x-axis, since 


(a.,, - a.) 
68 = — 1 1 - 1 


r. 


and 


So we have 


S'lT = r. (cos a i - cos a i+1 ). 
68 78W' = 1/c 2 , 


and we already know, of course, that 


S'/W' = 1/c 2 . 


Next, notice that the projection, 5UC of the (/+ l)-th arc on the y-axis equals 


80€ ' = r. (sin a 1+1 - sin c^) = c 2 [~ 


r 1 - 

- 1 1 

L T 

L i+i 

r i-i j 


The height at the tip of the i-th arc then is 

2X1 


1 


1 


= A- + - 

1 Lr i r i-l r 0 


"| for i = 1,2 ... n -1, 


since 


' = r Q sin Gj = c z /r v 


Also 


D6' 


[—- 1 ! 

Lr n-1 r 0 J 


+ r 


n-l 


Three, Four and Five Arcs 

The optimum solution for three arcs gives 


"X 


9 





36 = 1.399926... 


6 = 1.456879... » 0.9274780*(tt/2) 
f= 1.929128... 1.228121 *(n/2), 

and shows us that the ellipse is not optimal after all. For four arcs we find 

36 = 1.462089... 

6 = 1.448212... 25 0.9219604*(tt/2) 
y = 1.987501... 2S 1.265282*(w/2), 

and for five arcs we get 


36 = 1.500993 


• » • 


8 = 1.443930... » 0.9192345*(tt/2) 

= 2.024437... 25 1.288796>/2). 

We see that 8 is dropping more and more slowly, while If is growing, 
For five arcs, the parameters for the unsealed curve (r. = 1), are as follows 


r x = 0.500258... 

a x = 0.078707 

r 2 = 0.334497... 

a 2 = 0.237276 

r 3 = 0.253163... 

a 3 = 0.483045 

r 4 = 0.207337... 

a 4 = 0.847067 

r 5 = 0.189705... 

a 6 = 1.57079.. 


Perhaps we can guess the true minimum energy curve from the numerical data obtained so far. The 
parameters seem to roughly fit into a pattern like 


y 

These solutions are shown in Figure 5. 
as is 36. 


r - -—and n -TLllUli 

r i (/+1) and “i “ 2 «(«+!)' 


In this case 


(®i+i ~ «i) = w 


/+ 1 


n (n+1) 


so the arc lengths are 


&y r i( a i+i a i 


«(« + !) ' 


ft 

That is, the arcs all have the same length, and curvature increases linearly along the curve. 
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Now wc could sum the series for arc length and energy (easy) and for height and width (hard). Then we 
would discover that g = S ' W' decreases with n, and we could find its limit as the number of arcs tends to 
infinity. Instead, we proceed directly to the curve obtained in the limiting process. 


The Cornu Spiral 

The curve which has curvature varying linearly with are-length is called the Cornu Spiral (or Ruler’s Spiral) 
[pg. 190, 2; pg. 190, 14]. It can be defined using the two Fresnel integrals [pg. 820, 10; pg. 930, 11; pg. 300, 

12 ]. 


C(s) = f S cos (t! 2 ) dt, 
ao l 


If we let 


S(s) = J 5 sin (yf z ) dl. 
x - C(s) and y - S(s), 


we obtain a curve starting at the origin and curling upwards in the first quadrant. We note that 


x = cos (y.v 2 ) and y = sin (-jj*) 


7L.Z 


x — - wssin(?s z ) and y = 7 rscos(y/) 
This verifies that s is the arc-length along this curve, since 

x? + y z = 1, 

and that the curvature varies linearly with arc length, since 


7L.Z 


♦ • • • ♦ ♦ 


yx-xy = ws. 


The part of the spiral of interest to us here extends to the right up to the point where the curve becomes 
vertical, that is, x— 0. This is the point where s= 1 and 

jc = C(l) = 0.7798934... and y = S(l) = 0.4382591... 


The energy in this portion of the spiral is just 


g 


'= n k z ds- f {Tts) z ds = ir 2 /3. 

Jq J o 


We now build a smooth curve connecting the two points (-1,0) and ( + 1,0) by scaling, rotating, and shifting 
this tendril as shown in Figure 6. In the right hand quadrant we use 

Y — i _ SiA 3nf i v — 
x - 1 ^ and y - 

(The rest of the curve is obtained by reflection about the y-axis.) 


A" ■ ■ 



Wc then find that 


% = C(1)/,S’(1) = 1.779525... 

S = (tt 2 /3)A(1) = 1.441814... ^ 0.9178877*(tt/2), 

S = 1/.SY1) = 2.281755... ~ 1.452610*(w/2). 

The curve constructed out of a portion of die Cornu Spiral only has 91.78% of the energy of die semicircle 
and is thus die best curve so far. 

But, can we do better still? 


Six Arcs and More 


Unfortunately, the Cornu Spiral is not optimal either, as one sees by considering the best six-arc solution for 
which 

g = 1.441508... » 0.9176931 *(tt/2). 

For eight arcs, 
for sixteen arcs, 
for thirty-two arcs, 
and for sixty-four 

These solutions are shown in Figure 7. 

It seems that the total energy is approaching some limit, near 91.39% of that in the semicircle. Some of these 
results are summarized in Tables I and II. 


TABLE I 


n 

g/(7r/2) 

if/(7T/2) 

% 

r min 

1 

1.0 

1.0 

1.0 

1.0 

2 

.9415383 

1.160608 

1.284916 

.9245284 

3 

.9274780 

1.228121 

1.399925 

.8950140 

4 

.9219604 

1.265282 

1.462089 

.8794932 

5 

.9192345 

1.288796 

1.500993 

.8700129 

6 

.9176896 

1.305012 

1.527622 

.8636634 

7 

.9167300 

1.316870 

1.546987 

.8591353 

8 

.9160932 

1.325918 

1.561701 

.8557555 


g g ^ 0.916097*(7r/2), 
S 16 » 0.914532*(w/2), 
g 32 ^ 0.914285 *(tt/ 2), 
g 64 SS 0.913953*(m/2). 




TABLE!! 


n 

S/(tt/2) 

y/(w/2) 

% 

r mi n 

16 

.9144692 

1.358673 

1.614492 

.8442353 

32 

.9140406 

1.375704 

1.641625 

.8389106 

64 

.9139305 

1.384403 

1.655392 

.836516 

128 

.9139025 

1.388789 

1.662308 

.835461 

256 

.9138955 

1.390998 

1.665786 

.834998 

512 

.9138938 

1.392110 

1.667534 

.83473 

1024 

.9138934 

1.392671 

1.668417 

.83467 

2048 

.9138932 

1.392984 

1.668908 

.83463 


We can get better and better approximations to the optimum curve, provided we also carry out computations 
with more and more significant figures. Note, by the way, that while g varies little once n is reasonably large, 
y and % continue to show appreciable changes. This is a reflection of the fact that some distortions of the 
optimum curve produce only small changes in die total energy. 


Some Observations About the Optimum Curve 


The multi-arc approximation tends to the optimum curve in the limit as n tends to infinity. So we can learn 
some properties of the optimum curve from what we have so far. First of all, 

from r i r i-i s in a i = ° 2 

we get c 2 k 2 = cos 

where k is the curvature and <// is the angle which the curve makes with the x-axis. The constant c only affects 
die size of the curve, not its shape. We cannot determine it at diis point. 


From 


we get 


„ cos a,,, _cos a, 

r 2 -^- 1 


a i+l 


a i 


c z , 


„ sin [(a.,. —a .)/2] , 

'i (a. + 1 —a.)/2 sm K a i + l + ~ C ’ 


which in the limit again leads to 
We also obtain 


d& 

dx 


= cos 1 ^. 
= 1/c 2 , 


from 




Now 


g = f k z ds= Jk 2 -]/ 1 + (dy/dx) 2 dx. 






so 


k 2 -y 1 + (dyidx) 2 = 1/c 2 . 


Further, since, 
we again obtain, 

Each of these approaches leads us to the same 


* = “ *■ 

c 2 k 2 = cos y. 

simple differentia] equation for the curve. 


Finally, from 



we get in the limit 



That is, die curvature varies linearly along the axis of symmetry of die optimal curve. Substituting for k we 
also derive 


cos i p = (y/2c) 2 . 

Note that since the optimal curve bends downwards, its second derivative is negative. This is why, by the 
usual sign conventions, curvature too is negative. Thus we will use the equation 

— C K = V COS Ip 

between y = + 7 r /2 at die left end and y = — 77/2 at the right end of the curve. 




Differential Equation for the Curve 


The curvature is die rate of turning as one goes along the curve, diat is, 


K 


d 

' I i Jhii i i 

- ds ’ 


so 


and consequently, 


C^ = /cosy, 

f 

- C J7^J 


The substitution cos y = £ leads to a denominator of the form 


/id-iHi+i), 


while the substitution cos y = leads to a denominator of the form 

V l - £ 4 - 

In each case we are dealing with the integral of a rational function of £ and the square root of a cubic or 
quartic polynomial in £. This means that the answer can be expressed as an elliptic integral [pg. 16, 9; pg. 833, 
10; pg. 904,11; pg. 589,12], 
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Wc can actually just look up the result directly [using pg. 154,11; also see Appendix] and find the solution 

' / 

, 5 = V~2 c F(cos-Vcos?, 1 // 2 ), 

* ,* 

where V is the incomplete elliptic integral of the first kind. The constant of integration has been chosen so 
that ip~0 corresponds to s-0. The result came out positive because the integration goes from ^ = 0 to 
negative values of > p. For the half of the optimal curve in the negative quadrant, a minus sign must be 
attached. Some readers may notice the similarity of the solution found above, to dre equation for a pendulum 
[pg. 28, 9], swinging from -m/2 to + 7 t /2 (where ip is die angle from the vertical, while s corresponds to the 

time). 

We now have the equation of dre curve sought after in Whcwell form [pg. 4, 14], namely arc length as a 
function of tangential angle. We can immediately also rewrite it in Ccsaro form [pg. 4, 14], namely as a 
relation between arc length and curvature: 

s = VI c F( cos -1 (- c k), 1/ VI). 

Both of the forms given above are intrinsic equations for die curve [pg. 40,2], 

We can easily compute the length of the curve from an initial point at the top, where ip = 0, to the point on die 
jr-axis, where ip - — m/2 , since 


F{m/2, 1/ VI) = K{\/V2) 

where K is die complete elliptic integral of the first kind. Now [pg. 909,11], 


K( 17/2) = A"(sin (tt/4)) = VI f* 



where T is the gamma function [pg. 821,10; pg. 933,11; pg. 255,12], There is an infinite product [pg. 938,11] 
for T(l/4) 4 


f(l/4) 4 = 16 m 2 


oo 


n 

ir - 1 


<4k-]) z \(4k+l) z -l] 
(4k+ 1) 2 [(4&-1) 2 -1] 


which gives us the numerical value 


The arc length is finally, 


r(l/4) = 3.6256099082... 



crnMii 

2 V 2m 


1. Note that some authors [pg. 833,10] list the arguments of the incomplete elliptic functions in the reverse 
order of that shown here. 



Cartesian Form of the Solution 


For many purposes it is more convenient to express the solution as a relationship between the x and y 
coordinates. 


We note that 


- cos and Q = sin 


and remember that 


cf = SZSTf. 


By the chain-rule for differentiation tt = ~r / = - cV cos 4, 

ay as 1 as T 


and 


% ~ fc/^ds = ~ csin ^ ^ cos 


The latter equation is easy to integrate using the substitution z = cos ip. 


y-c 


I-77 = 


so that 


(y/2c) = V cos \p. 


where the constant of integration was chosen so that>’=0 when \p = -tt/2. The reader may also recall that in 
the previous section this result was found directly as the limit of height of the multi-arc approximation. 


The integral for x is a little harder, 


X = C J V cos ip dip. 


This can be expressed as the difference of two incomplete elliptic integrals [using pg. 156, 11; also see 
Appendix], 


1 


* = V~2 c[ 2 Ffcos'V cos \p, 1//2) - F (cos' V cos ip, 1//2)], 
where the constant of integration was chosen so that x = 0 when ip = 0. Finally then. 


x = Vic [2 £(cos' 1 (y2c), 1/ VI) - F(cos“V2c), V VI) ]. 


An alternate way to obtain the same result is to note that x and y arc related by the differential equation 


dy _ V 1 - (y/2c)' 

dx ~ (yntf 


1. Note that E denotes two things in this paper. When the letter appears with two argument it signifies the 
incomplete elliptic integral of the second kind. 

2. Note again that some authors [pg. 833,10] list the arguments of the incomplete elliptic functions in the 
reverse order of that shown here. 
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sincc 


tan ip 


1 -cos 2 ^ 
cos ip 


In any case, for j; = 0, 


x = n c [ 2 /:(1 /V"2) — K( 1/ V”2) ], 


where A" and fare the complete elliptic integrals of the first and second kind respectively. 
Using Legendre’s identity [pg. 25, 9; pg. 836,10; pg. 907.11; pg. 591f 12] 


where 


we get [pg. 25, 9] 


so that 


IXk)K(k') + E(k') K(k) - K(k)K(k') = tt/2, 

k z + k' z - 1, 

2 /-;(1//2) = [ w/2 + A 2 (l//2) ] / A(l//1), 
[ 2 £’(l//2) - A'(l//2) ] = (tt/2) / £(1/71). 


So the width of the curve is 


W= c 


3/2 


r(i/4>‘ 


and, if we want W = 1, we must have 


_ _ 0.8346268416... 

(2tt) 3/2 


The height of die curve then comes to 


% = 2 c = 1.669253683... 


(compare to the two-arc approximation). The minimum radius of curvature, the inverse of the maximum 
curvature, is 


_ — L— = c = 0.8346268416... 
(%/2 c z ) 


Thus a circle tangent to the curve at the top is also tangent to the x-axis. 


The arc length comes to 


Finally, from 


we get 


and so 


Note that 8 if 


t - 1 rtl/4) 4 

)mes to y - 2 {2in)Z , 

y = 2.188439615... « 1.393203929*(w/2). 

s=^ 

8 = 1/c 2 , 

g = -^^4 = 1.435540022... ~ 0.9138931623 *(tt/2). 
77 . The curve of least energy is shown in Figure 8. 




Extension of the Curve 


So far we have considered a finite segment of the optimum curve, extending from a point of zero curvature 
(where 4'= -v/2) through a point of maximum curvature (where i// = 0) to a second point of zero curvature 
(where tp = + u/ 2). Can the curve be extended beyond these points? 

It is clear that $ must remain in die range [-w/2, + 7 t/ 2] so tliat the square root of its cosine remains real. 
To continue the curve then, the sign of the curvature must change; we must choose the other sign for the 

square-root. The new segment we obtain has the same shape, of course, as die segment we have found 
already, just inverted. 


The segment we have used so far is just a piece of the infinite periodic wave shown in Figure 9. In Figure 10 
we see several curves which correspond to stationary values of the integral and pass through die specified 
points with the desired orientation. The one on the left is the one which corresponds to a global minimum of 

the energy. The curves containing n half-cycles have an energy n 2 as large as the one containing a single 
half-cycle. 

The curve of least energy passing through two given points with specified orientation is just a portion of the 
general curve, suitably translated, rotated, and scaled. This is illustrated in Figure 11. The rotated, translated 
and scaled curves form a four parameter family. 


Variational Approach 

We are trying to find die curve for which 


f rZ ^ s 


is minimal. This integral can also be written in the form 


8 = JV [1 + (yf] 1/z dx. 


or, since 


K 


JL 


n 


[i + 001 


'\2i3/2 > 


as 


g 


= f 


"\2 


in 


[1 + Oil 


'\2i5/2 


dx 


ITiis is of the form 


f T (x, y, y, y") dx. 


and the calculus of variation [pg. 119,15; pg. 190,16; pg. 198,17] teaches us that for a stationary value of the 
integral, 

*y ~ i*Y + &*? = *■ 

where 9 y , 5 », and SF" are the pardal derivatives of If with respect toy', y, and y" respectively: 
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‘ J y = o, 


5 /O") 2 

[1 + (y) z ] 1/z ' 


V [1 + 0') 2 ] 5/2 


Since = 0, 


*[*?’-*v] =0 - 


and integrating, we get 


V + & V 


where /I is an arbitrary constant. In the above we have closely followed the approach taken by Mehlum 
[pgs. 157&189,2; pg. 43, 6]. 


Also 


jL 9 2 r 10 

dx*y [1 + (y) z f n [1 + (y') z ] 1/z ’ 


so we get 


Now 


[l + 0') 2 1 5/2 [1 + (y'Y) 


X k") 2 .- 

/.a2i7/! 


dy' d£_ ,dy_ . „ 

-ir =az= y f y • 


// dy 

y =y aJ • 


In addition 


JL 1 d£_ 

dy 2 (y) ~ y dy’ ’ 


and 


<*/ [l + O ') 2 ] 5 7 2 


[i + or] 


'\2i7/2 


Using these results in the equation above, 


,"\2 


[l + o') ] 


l _J 

a2i5/2 “ ^ » 


,"\2 


and, integrating, we get 

where Z? is a second arbitrary constant. 


[i 7WF* ~ Ay ' + B ■ 


Returning for a moment to the integral of the square of curvature, we see that 
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J n z ds — J [Ay + B] dx — J A dy + J B 
— A O'j — j' 0 ) + B(x 1 — x Q ) 


dx 


for a curve which starts at (x Q , ;- Q ) and ends at (x : , y ). 


Now 


so 


(£l 


[1 + (/) ] 


' v 2i5/2 


k Z - y / 1 + (/) 2 , 


k Z ~\J 1 ■+■ (y ) 2 — Ay + B. 


Also. 


dy dy , dx .,. 
fx = dsljs =y/x ' 


so 


or 


since 


Now if 


then 


Remembering that 
we finally see that 
Letting 


we get 


k 2 -y/ 1 + (y/x) 2 = A (j’/x) + B , 


ic 2 = Ay + Bx, 
x 2 + y 2 = 1. 


tan \p 


dy 

dx' 


x = cos 4 1 and y = sin . 


K 


d± 
ds ’ 


d± 

ds 


± V /I sin + 5 cos . 


/I = 1/c 2 cos <p and B - — 1/c 2 sin <p 




ds 


± V sin (<f — <p). 


The scale of the curve is dependent on the parameter c, while the rotation in the xj-plane is dependent on the 
parameter y. 


Altogether we have a four parameter family of curves, since we can also choose an initial point and a direction 
for the curve. Conversely, we can find a single curve out of this family which passes through any two points 
with orientation specified at both points. 

By the way, if we let the line from the initial point (* 0 j' 0 ) to the final point (x^) have length r and 
direction 6, then 


JK l ds= (r/c 2 ) sin (0 - <p). 
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Summary 

Wc have found the simple equation 

± c k - V cos (rp - q >), 

and solved it to find the equation of die curve of least energy in Whewell form 

s = ± V~2 c F( cos’ 1 -/cos (ip — fp), 1/V2). 


We also developed a differential equation for the curve in Cartesian coordinates aligned with the axis of 


symmetry, 

d}> _ _ V 1 ~ Q/2c) 4 

dx ~ (y/2cf 

The solution of this equation, for given initial conditions, led to 

* 

x - -/2c [ 2 E (cos‘ 1 (y/2c), 1/V~2) - / r (cos' 1 (>72c), 1/-/2)]. 


Wc considered the curve of least energy connecting the point ( — 1,0) to the point ( + 1,0) with vertical initial 
and final orientations. This curve has minimum radius of curvature 

_ rn/4) 2 
C ” (2tt) 372 ’ 


rises to a height 
has arc length 

and energy 


% = 2c, 

1 10 / 4)1 

2 (2tt) 2 

c _ -12 7r ) 3 

T(l/4) 4 ’ 


or about 91.39% of that of the simple semicircle approximation. 

We have also given a method for finding approximations, consisting of circular arcs, to the curve of least 
energy. 

Note that the curve found here is extensible [7] in the sense that if the least energy curve with orientation a at 
A and orientation ft at B passes through the point C with orientation y, then the segments from A to C and 
from C to B are themselves least energy curves. (This is not true of the two-arc approximation [8]). As a 
result, such a curve can be computed by a simple, locally connected network. 

We have not shown how to find the particular member of the four parameter family of curves which passes 
through a specified pair of points with specified orientation. Presumably determining the axis of symmetry 
of the curve would be a helpful first step in this direction. In practical applications the multi-arc 
approximation method may be suitable in dealing with this problem. We have not discussed how one might 
compute the curve of least energy passing through three or more points. Here there is no constraint on the 
direction, but the curvature must be continuous. Nor have we touched upon the extension to curves and 
surfaces in three dimensions, a topic which Mchlum addresses [pg. 62,6]. 
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Appendix 

In this appendix we determine two integrals needed in the body of the paper. 



To evaluate the integral 



we substitute 


and obtain 


Next, substituting 


we get 


and since 



we finally get 


cos 


= /2 Ficos’ 1 V cos \p, l/V~2). 



To evaluate the integral 



V cos if dip, 



The integral thus can be split into two parts, the second of which we have already evaluated. 
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2 2 

we substitute as before r — 1 - £ 


and obtain 


and so 

since sin" 1 V 1 - cos \p = cos' 1 -/cos . 



/„-/.= 2/~2 /i (cos' 1 / cos /, 1/-/2), 


Finally then, 




-/2 [ 2 E (cos' 1 V cos /, 1//~2) - F(cos* 1 V cos 1/-/"2) ]. 


The corresponding definite integrals from ^ = 0 to / = w/2 can be expressed in terms of complete elliptic 
integrals of the first and second kind [pg. 91, 9], 
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Figure ]: A semicircle connects the points (-1,0) and (+1,0) and has vertical orientation at these points. The 
arc in the first quadrant has length w/2and energy n/2. 































Figure 4: A curve constructed out of n circular arcs. Each arc contributes to the total arc length and the to 
energy. 

















Figure 6: A curve constructed from part of the Cornu Spiral. Here the curvature varies linearly with arc 
length along the curve. 































Figure 8: The curve of least energy. Here the curvature varies linearly with distance along the axis of 
symmetry. A circle tangent at the top is also tangent to the jt-axis. 


















etween — *n/2 and 4*7r 










Figure 10: Curves corresponding to stationary 
cycle, (B) 3/2 cycles, (C) 5/2 cycles, 












Figure 1]: Optimum curve passing through two given points with given orientation. It is a member of a four 
parameter family of curves obtained by translating, rotating, and scaling the particular curve 
passing through the points (— 1,0) and ( + 1,0) with vertical orientation. 








